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' Abstract 



We study unitary random matrix ensembles of the form Z^^jyl detMp"e~^'^''^'*^''c?M, 
where a > —1/2 and V is such that the Hmiting mean eigenvalue density for n, N oo 
and n/N 1 vanishes quadratically at the origin. In order to compute the double scaling lim- 
' its of the eigenvalue correlation kernel near the origin, we use the Deift/Zhou steepest descent 

I method applied to the Riemann-Hilbert problem for orthogonal polynomials on the real line 

■ with respect to the weight \x\'^"'e~^^^^\ Here the main focus is on the construction of a local 
parametrix near the origin with i/j-functions associated with a special solution Qa of the Painleve 

■ II equation q" = sq + 2q^ — a. We show that has no real poles for a > —1/2, by proving the 
lO ! solvability of the corresponding Riemann-Hilbert problem. We also show that the asymptotics 

' of the recurrence coefficients of the orthogonal polynomials can be expressed in terms of qa in 

^ I the double scaling limit. 

"5 '. 1 Introduction and statement of results 
1.1 Unitary random matrix ensembles 

For n G N, > 0, and a > —1/2, we consider the unitary random matrix ensemble 

Z-jv|detM|2-e-^Try(Af)^^^ ^^^^^ 

on the space of n x n Hermitian matrices M, where y : M ^ M is a real analytic function satisfying 



lim —^^f)-^ = +oo. (1.2) 

x^±oo log(x^ + Ij 

Because of (|1.2() and a > —1/2, the integral 

Z„,^ = y"|detM|2"e-^Ttv(M)^^ ^^3) 

converges and the matrix ensemble ()1.1() is well-defined. It is well known, see for example jHHSOl) that 
the eigenvalues of M are distributed according to a determinantal point process with a correlation 
kernel given by 

n-l 

KnM^,y) = |xre-f^(")|yre-f^(^)^Pfc,,v(x)pfc,iv(y), (1.4) 

k=0 
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where Pk,N = i^k nx^ + ■ ■ ■ j i^k,N > 0) denotes the k-th degree orthonormal polynomial with respect 
to the weight |x|'2"e-^^(^) on M. 

Scaling limits of the kernel H1.4() as n, ^ oo, n/N 1, show a remarkable universal behavior 
which is determined to a large extent by the limiting mean density of eigenvalues 

ipvix) = lim -Kn,nix,x). (1.5) 

Indeed, for the case a = 0, Bleher and Its U| (for quartic V) and Deift et al. fT3^ (for general real 
analytic V) showed that the sine kernel is universal in the bulk of the spectrum, i.e., 

1 / u V \ simriu — v) 

hm — — — -Kn,n xo H — — -, xo + 



n^oo r^^^/(xo) ' V '>^'^lJv{xo)^ ml^vixo) J ■7t{u — v) 

whenever ipvixo) > 0. In addition, the Airy kernel appears generically at endpoints of the spectrum. 
If xq is a right endpoint and ipv{x) ~ {xq — x)^^"^ as re ^ xq—, then there exists a constant c > 
such that 

1 / u V \ Ai(u)Ai'(v) - Ai'(u)Ai(v) 



n— »oo cn^''^ V cn^/'^ cn^/'^/ u — v 



where Ai denotes the Airy function, see also |E 

The extra factor |detMp" in introduces singular behavior at if a 7^ 0. The pointwise 
limit (HHl) does not hold if V'y(O) > 0, since Kn^ni^, 0) = if a > and Kn^n{0, 0) = +00 if a < 0, 
due to the factor in ()1.4|) . However H1.5() continues to hold for x ^ and also in the sense 

of weak* convergence of probability measures 

1 * 

—Knnix,x)dx^ipY(x)dx, as n — > cxD. 
n 

Therefore we can still call i/jy tlie limiting mean density of eigenvalues. Observe that ipy does not 
depend on a. 

However, at a microscopic level the introduction of the factor | det Mp'^ changes the eigenvalue 
correlations near the origin. Indeed, for the case of a non-critical V for which ipvi^) > 0, it was 
shown in (29j that 

lim r^K, 



TrVuVv ^^77 , 1-6 

2{u — V) 

where Jj^ denotes the usual Bessel function of order f. 

We notice that universality results for orthogonal and symplectic ensembles of random matrices 
have been obtained only very recently, see [01 E3 ^3 • 

1.2 Multi-critical case 

It is the goal of this paper to study in a critical case where tpv vanishes quadratically at 0, 



I.e., 



V'y(O) = V'y(O) = 0, and Vy(0) > 0. (1.7) 
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The behavior (jl.Tf) is among the possible singular behaviors that were classified in ^2]- The 
classification depends on the characterization of the measure ipv{x)dx as the unique minimizer of 
the logarithmic energy 

Iv{^A =11 \og—^d^l{x)d^l{y) + / v{x)d^Ji{x) (1.8) 



\x - y\ 

among all probability measures on M. The corresponding Euler-Lagrange variational conditions 
give that for some constant £ G M, 



log \x — y\ipviy)dy — V{x) + £ = 0, for x G supp(V'y), (1-9) 

2 j log \x - y\i^v{y)dy - V{x) + ^ < 0, for x G R. (1.10) 

In addition one has that is supported on a finite union of disjoint intervals, and 

i^vix) = l^Qvix), (1.11) 

where Qy is a real analytic function, and Qy denotes its negative part. Note that the endpoints 
of the support correspond to zeros of Qy with odd multiplicity. 
The possible singular behaviors cirG clS follows, see |12| I26j. 

Singular case I ipv vanishes at an endpoint to higher order than a square root. This corresponds 
to a zero of the function Qy in (|1.1H) of odd multiplicity 4fc + 1 with k > 1. (The multipicity 
4/c + 3 cannot occur in these matrix models.) 

Singular case II ipv vanishes at an interior point of supp(V'y), which corresponds to a zero of 
Qv in the interior of the support. The multiplicity of such a zero is necessarily a multiple of 
4. 

Singular case III Equality holds in the variational inequality (jLlOj) for some x G M \ supp(^/'y). 

In each of the above cases, V is called singular, otherwise regular. The above conditions correspond 
to a singular endpoint, a singular interior point, and a singular exterior point, respectively. 

In each of the singular cases one expects a family of possible limiting kernels in a double scaling 
limit as n, ^ oo and n/N ^ 1 at some critical rate jSj. As said before we consider the case (|1.7|) 
which corresponds to the singular case II with = 1 at the singular point x = 0. For technical 
reasons we assume that there are no other singular points besides 0. Setting t = n/N, and letting 
n, — > cxD such that t ^ 1, we have that the parameter t describes the transition from the case 
where il'viO) > (for t > 1) through the multi-critical case (i = 1) to the case where lies in a gap 
between two intervals of the spectrum (t < 1). The appropriate double scaling limit will be such 
that the limit liuin^N^oo n"^^^ {t — 1) exists. 

The double scaling limit for a = was considered in jJl IE] for certain special cases, and in 
[7] in general. The limiting kernel is built out of -^-functions associated with the Hastings- McLeod 
solution j2Qj of the Painleve II equation q" = sq + 2g^ . 

For general a > —1/2, we are led to the general Painleve II equation 

q" = sq + 2q^ - a. (1.12) 

The Painleve II equation for general a has been suggested by the physical papers The 
limiting kernels in the double scaling limit are associated with a special distinguished solution of 
(|1.12() . which we describe first. We assume from now on that a ^ 0. 
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1.3 General Painleve II equation 

Balancing sq and a in the differential equation ()1.12() . we find that there exist solutions such that 



a 



q(s) ~ — , as s — > +00, (1-13) 
s 

and balancing sq and 2q'^, we see that there also exist solutions of ()1.12|) such that 

qis) ~ \J^, as s^ -00. (1.14) 

There is exactly one solution of (|1.12j) that satisfies both (|1.13|) and l|1.14j) , see |2H 1221 , and we 
denote it by qa- This is the special solution that we need. It corresponds to the choice of Stokes 
multipliers 

si = e-™, S2 = 0, S3 = -e™, 

see Section 2 below. We call qa the Hastings-McLeod solution of the general Painleve II equation 
(|1.12() . since it seems to be the natural analogue of the Hastings-McLeod solution for a = 0. 

The Hastings-McLeod solution is meromorphic in s (as are all solutions of (|1.12j) ) with an 
infinite number of poles. We need that it has no poles on the real line. From the asymptotic 
behavior (|l.ll-{j) and (|1.14jl we know that there are no real poles for |s| large enough, but that does 
not exclude the possibility of a finite number of real poles. While there is a substantial literature on 
Painleve equations and Painleve transcendents, we have not been able to find the following result. 

Theorem 1.1 Let q^ be the Hastings-McLeod solution of the general Painleve II equation H1.12() 
with a > —1/2. Then qa is a meromorphic function with no poles on the real line. 

1.4 Main result 

To describe our main result, we recall the notion of ^/^-functions associated with the Painleve II 
equation, see jl6j . The Painleve II equation H1.12() is the compatibility condition for the following 
system of linear differential equations for ^ = ^a{C: s). 

_ = AW, - = B*, (1.15) 

where 

4_M<'-^(« + V) 4Cq + 2ir + a/C\ r - q\ n Ifi^ 

That is, H1.15() has a solution where q = q{s) and r = r(s) depend on s but not on C, if and only if 
q satisfies Painleve II and r = q' . 
Given s, q and r, the solutions of 

^r^.(c)^_4J.(c)^ 



dc \MC)J V^2(c) 
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are analytic with branch point at C = 0. For a > —1/2 and s G M, we take q = qa{s) and r = q'^[s 

'^a,i(C;s) 



where q^ is the Hastings-McLeod solution of the Painleve II equation, and we define 
as the unique solution of (|1.17|) with asymptotics 



'^a,2(C;s) 



„S,..,o(|»..(C; J)). 

uniformly as C — > cxd in the sector e < arg C < vr — e for any e > 0. Note that this is well-defined 
for every s S M because of Theorem ll.il 

The functions ^a,i and ^a,2 extend to analytic functions on C\ (— ioo,0], which we also denote 
by ^a,i and ^a,2, see also Remark 12 . 331 below. Their values on the real line appear in the limiting 
kernel. The following is the main result of this paper. 

Theorem 1.2 Let V be real analytic on M such that (|1.2|) holds. Suppose that ipy vanishes quadrat- 
ically in the origin, i.e., ipv{^) = V'y(O) = 0) o,nd ipy{^) > 0, and that there are no other singular 
points besides 0. Let n, N ^ oo such that 

lim n^/^{n/N -1) = L eM. 

n.N^oo 



exists. Define constants 



(1.19) 

and 

s = 27T'/^L[^lj'{.i0)]-'/'wsA0), (1.20) 
where wsy is the equilibrium density of the support of ipv (see B,emark \l .'A\ below). Then 

hm -l^Kr.^^(-!^,-^)=K'^"t^-{u,v;s), (1.21) 



n,N-^co cn^/^ ' \cn^/^' cn^/^ 
uniformly for u,v in compact subsets o/M\{0}, where 

K^rtt,a^^ ^. ^) ^ _gi™[sgn(u)+sgn(^)] ^a,l(^; g)^a,2(^; s) - ^a,l{v\ s)$a,2(«; s) ^2) 

' ' 2Tri{u — v) 

Remark 1.3 The equilibrium measure of Sy = supp{ipv) is the unique probability measure uJSy 
on Sy that minimizes the logarithmic energy 

Hl^) =11 log-^ — ^- — -dfi{x)dfi{y) 
JJ \x-y\ 

among all probability measures on Sy- Since Sy consists of a finite union of intervals, and since 
is an interior point of one of these intervals, has a density wsy with respect to Lebesgue 
measure, and wsyi^) > 0- This number is used in H1.2U() . 
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Remark 1.4 One can refine the calculations of Section|l]to obtain the following stronger result: 
' ^n.^(-4H>-47^)=^^"*'"(«>-;^) + of^S^V (1-23) 



uniformly for u, v in bounded subsets of M \ {0}. 

Remark 1.5 It is not immediate from the expression (|r^ that i^crii,^ -g rpj^-g pj-Qpgrty 

follows from the symmetry 



gi^sgnH^^^^^^. ^ e^™"S'^(")$„,i(u; s), forn G M \ {0}, 

which leads to the "real formula" 

K''"''"{u,v;s) = — -Im ( el™fe'iW-'^g°(«))$^ ^(u; s)$„ ^{v; s)) , 

7r[u — V) \ ' ' / 

see R emark 12.111 below . 

Remark 1.6 For a = 0, the theorem is proven in 0. The proof for the general case follows along 
similar lines, but we need the information about the existence of qa{s) for real s, as guaranteed by 
Theorem 11.11 

1.5 Recurrence coefficients for orthogonal polynomials 

In order to prove Theorem 11.21 we will study the Riemann-Hilbert problem for orthogonal poly- 
nomials with respect to the weight \x\'^°'e~^^^^\ This analysis leads to asymptotics for the kernel 
Kn^N, but also provides the ingredients to derive asymptotics for the orthogonal polynomials and 
for the coefficients in the recurrence relation that is satisfied by them. 

To state these results we introduce measures vt in the following way, see also and Section 3.2. 
Take Sq > sufficiently small and let be the minimizer of ly/^^v) (see (|1.8|) for the definition of 
/y) among all measures u = v'^ — where are nonnegative measures on M such that i/(M) = 1 
and supp(z^^) C [—5q, 5q]. We use V'* to denote the density of f^. 

We restrict ourselves to the one-interval case without singular points except for 0. Then 
supp(V'y) = [fl) b] and supp{^pt) = [a*, h] for t close to 1, where at and bt are real analytic functions 
oft. 

We write TTn^N for the monic orthogonal polynomial of degree n with respect to the weight 
^^^2a^-NV{x) ^ Those polynomials satisfy a three-term recurrence relation 

T^n+l,N = - hn,N)T^n,N - a^,N'^n-l,N, (1-24) 

with recurrence coefficients an,N and hn,N- In the large n expansion of an,N and bn,N-, we observe 
oscillations in the C'(n~^/^)-term. The amplitude of the oscillations is proportional to Qais), while 
in general the frequency of the oscillations slowly varies with t = n/N . + 

Theorem 1.7 Let the conditions of Theorem, 11.21 he satisfied and assume that supp(?/'v') = V^^b] 
consists of one single interval. Consider the three-term recurrence relation H1.24() for the monic 
orthogonal polynomials T^k^N with respect to the weight |xp°e~^^^^-* . Then as n, N ^ oo such that 
n/N — 1 = 0{n^'^/^), we have 

4 2c 
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2 c 
where t = n/N , c is given by l\l.W\i . 

st,n = n^/^-MO), (1.27) 



6* = arcsin^-^, (1.28) 
b — a 

and 

0Jt= f\t{x)dx. (1.29) 



Remark 1.8 It was shown in [7] that ^^f(0)|j_-^ = wsy{^), which in the situation of Theorem 
11.71 imphes that (since Sy = [a, b] and ^t(O) is real analytic as a function of t near t = 1), 

V^t(O) = (t-l)— i= + 0((t-l)2), ast^l. 
vrv — ao 

Then it follows from (|1.27jl that St^n = n^^^it — 1) ^^/^-^ + 0{n^'^^'^) and we could in fact replace 
st,n in (fT^ and (fOH)) by 

1 



c\/—ab 



We prefer to use st^n since it appears more naturally from our analysis. 

Remark 1.9 In j5|, Bleher and Its derived ()1.25p in the case where a = and where ^ is a critical 
even quartic polynomial. They also computed the C'(n~^/^)-term in the large n expansion for an,N- 
For even V we have that a = —b, 9 = 0, ujt = 1/2 and thus cos(27rna;t + 2a6) = (—1)", so that 
(|1.25|) reduces to 



b _ g^Kn)(-l)\ ~i/3 + o(n-2/3), 



On M 

' 2 2c 

which is in agreement with the result of [Sj. Also for even V the recurrence coefficient bn,N vanishes 
which is in agreement with (|1.26j) . 

Remark 1.10 In 3^ an ansatz was made about the recurrence coefficients associated with a general 
(not necessarily even) critical quartic polynomial V in the case a = 0. For fixed large A^, the ansatz 
agrees with (|1.25|) and (|1.26jl up to an A^-dependent phase shift in the trigonometric functions. 



1.6 Outline of the rest of the paper 

In Section 2, we comment on the Riemann-Hilbert problem associated with the Painleve II equation. 
We also prove the existence of a solution to this RH problem for real values of the parameter s, 
and this existence provides the proof of Theorem ll.il In Section 3, we state the RH problem for 
orthogonal polynomials and apply the Deift/Zhou steepest descent method. Our main focus will 
be the construction of a local parametrix near the origin. For this construction, we will use the RH 
problem from Section 2. In Section 4 and Section 5 finally, we use the results obtained in Section 
3 to prove Theorem 11.21 and Theorem 11.71 
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Ta r4 
Figure 1: The contour S consisting of four straight rays oriented to infinity. 

2 RH problem for Painleve II and the proof of Theorem 11.11 

As before, we assume a > —1/2. 

2.1 Statement of RH problem 

Let S = \Jj Tj be the contour consisting of four straight rays oriented to infinity, 

vr 57r 5tt vr 

ri:argC = -, r2 : argC = — , r3:argC = — —, r4:argC = --. 

DO DO 

The contour S divides the complex plane into four regions Si, . . . , S/^ as shown in Figure ^ For 
a > —1/2 and s e C, we seek a 2 x 2 matrix valued function ^'q,(C;s) = ^'a(C) (we suppress 
notation of s for brevity) satisfying the following. 

RH problem for ^'^r 

(a) is analytic in C \ S. 

(b) satisfies the following jump relations on S \ {0}, 

*„,+ (C) = ^a,-(C)(^-L forCGTi, (2.1) 

vI'«,+ (C)=^„,-(C)(_gL i), forCGTs, (2.2) 
*«,+ (C)=^„,-(C)(^0 1 J, forCGTs, (2.3) 

^'„,+ (C) = ^a,-(C)(^0 \ j, forCGr4. (2.4) 

(c) has the following behavior at infinity, 

^'„(C) = (/ + 0(l/C))e-'(^'^''"'^^"^ asC^oo. (2.5) 
Here o"3 = ( o -i ) denotes the third Pauli matrix. 
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(d) has the following behavior near the origin. If a < 0, 



and if a > 0, 



o 



o 



o 



a 


ICI \ 




ICIV ' 


ici- 




ici- 




icr 


icr" 


icr 


icr" 


ici- 




ici- 





as C ^ 0, 



(2.6) 



(2.7) 



, as C ^ 0, C e 5i U ^3, 
as C ^ 0, C G 5*2, 
as C — > 0, C £ 5'4. 

Note that ^'q, depends on s only through the asymptotic condition 1)2. 5|) . 

Remark 2.1 This RH problem is a generalization of the RH problem for the case where a = 0, 
used in [im. 

Remark 2.2 By standard arguments based on Liouville's theorem, see e.g. jHlIlTj, it can be verified 
that the solution of this RH problem, if it exists, is unique. Here it is important that a > —1/2. 

In the following we need more information on the behavior of solutions of the RH problem near 
0. To this end, we make use of the following proposition, cf. P2?- We use Gj to denote the jump 
matrix of on Tj as given by (|TT |) - ((n)) . 

Proposition 2.3 Let ^ satisfy conditions (a), (b), and (d) of the RH problem for 
(1) If a — ^ ^ No, then there exists an analytic matrix valued function E and constant matrices 



Aj such that 



*(C) = EiC) 







Aj, forCeSj, 





where the branch cut of is chosen along T4 . The matrices Aj satisfy 
Aj+i = AjGj, for j = 1,2,3, 

and 

'0 



(2i 



(2.^ 



A. 



P 



-1 



2 cos(7ra) 



for some p € C \ {0}. 



(2.10) 



(2) If a — ^ ^ No, then there is logarithmic behavior of at the origin. There exists an analytic 
matrix valued function E and constant matrices Aj such that 



*(C)=i?(C)( i^Cj"^^ forCeSj, 



(2.11) 
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where again the branch cuts of and ln( are chosen along The matrices Aj satisfy 

Aj+i = AjGj, for j = 1,2,3, (2.12) 
and for some p £ C, 

if a — ^ is even, 

(2.13) 

if a — ^ is odd. 

Proof. (1) Define E by equation (|2.8|) with matrices Aj satisfying 1)2. 9|) and ()2.1U|) . Then E is 
analytic across Fi, and because of 1)2. 9() . For C ^ ^4, there is a jump 

E+ic) = E.io l^^'J ^4G4Ar'(^o" A)^- ^^-^^^ 

Using = e^'^*"(^" and the exphcit expressions for the matrices Gj and Aj, we get from (|2.14() 
that E is analytic across r4 as well. 

What remains to be shown is that the possible isolated singularity of E at the origin is removable. 
If Q < it follows from (|2IEI) and (|2IHI) that 

^(C) = o(j^j2l }), asC-0, 

so that (since 2q > — 1) the isolated singularity at the origin is indeed removable. If a > we have 
in sector S2 by ^1^, ((THl), and (I^TTTl) that 

as C — > 0, C G S2, where * denotes an unimportant constant. Hence the singularity at the origin is 
not a pole. Moreover, from ()2.7() and (|2.8() it is also easy to check that E does not have an essential 
singularity at the origin either. Therefore the singularity is removable for the case a > as well, 
and the proof of part (1) is complete. 

(2) The proof of part (2) is similar. □ 

Remark 2.4 The matrix A2 in Proposition 12 .31 is called the connection matrix, cf. jlbl I19j . In all 

cases we have det ^2 = 1 and the (1, l)-entry of A2 is zero. 

2.2 Solvability of the RH Problem for ^„ 

We are going to prove that the RH problem for is solvable for every s G M. We do that, as 
in jl3l 1191 by showing that every solution of the homogeneous RH problem is identically zero. 
Such a result is known as a vanishing lemma |18[ I19j. 

We briefly indicate why the vanishing lemma is enough to establish the solvability of the RH 
problem for The RH problem is equivalent to a singular integral equation on the contour S. 
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The singular integral equation can be stated in operator theoretic terms, and the operator is a 
Fredholm operator of zero index. The vanishing lemma yields that the kernel is trivial, and so the 
operator is onto which implies that the singular integral equation is solvable, and therefore the RH 
problem is solvable. For more details and other examples of this procedure see |131 1191 136j and |231 
Appendix A]. 

Proposition 2.5 (vanishing lemma) Let a > —1/2 and s G M. Suppose that ^ satisfies condi- 
tions (a), (b), and (d) of the RH problem for ^ ^ with the following asymptotic condition (instead 
of condition (c) ) 



O{l/0, 



as C ^ oo. 



(2.15) 



Then = 0. 



Proof. As before, we use Gj to denote the jump matrix of Tj, given by (|2.H) - (|2.4|) . Introduce an 
auxiliary matrix valued function H with jumps only on M, as follows. 



H{0 



for C G ^2 U 54, 
for C G 5i n C+ 



^(C)Gie^(3C'+<)'^3, 
$(C)G2 ^e*(t^'+<)'"3^ for (eSsD C+, 

for C G 53nC_, 



for C G 5i n I 



Then H satisfies the following RH problem. 

RH problem for H: 

(a) : C \ M — > C^^^ is analytic and satisfies the following jump relations on M \ {0}, 



-e 



^ je*(3C'+<)-3^ for Cg (-cx),0), 
e^(tf'+<)-3^ forCG(0,( 



, oo 



(b) H{C) = 0{1/C), asC-oo. 

(c) H has the following behavior near the origin: If a < 0, 
and if a > 0, 



as C — > 0, 



H{C) = 



O 



o 



icr 


ici-" 


icr 


icr" 


ici-" 


icr 


ici-" 


icr 



as C ^ 0, Im C > 0, 
as C ^ 0, Im C < 0. 



(2.16) 



(2.17) 
(2.18) 



(2.19) 



(2.20) 
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The jumps in (a) follow from straightforward calculation. The vanishing behavior (b) of H at 
infinity (in all sectors) follows from the triangular shape of the jump matrices Gj, see H'2.ip - ()2.4[) . 
For example, for C G 5"! n C+ we have Rei(|C^ + sC) < so that by p?T5|l and (j2l6|l 



HiC) = 0(1/0 



1 







g-™g2i{|C^+sC) I 



o{i/C), 



as C — > oo. 



The behavior near the origin in (c) follows from ProDosition l2.31 This is immediate for (|2.19|) . while 
for a > 0, a - i No, we have by (|THI), (HSJ, (I^Tini) . and (ITTKl) . 



E{C)C---^ A2 = EiOC''' \ * ) , if Im C > 0, 

\ * * / 



i?(C)C-""^^4 = i?(C)C-""^ ( * ° ) , if ImC < 0, 



which yields H2.20() in case a — | No, since E is analytic. Using (|2.13() instead of 1)2. lOj) . we will 
see that the same argument works if a — | S No. 

Next we define, cf. [13113120!, 



MiC) = HiOHiCr, forCGC\R, 



(2.21) 



where H* denotes the Hermitian conjugate of H. From condition (c) of the RH problem for H it 
follows that M has the following behavior near the origin 



M(C) = <^ 



^l|^|2a |^|2a)' as ( ^ 0, m casc « < 0, 



O 



1 1 



as C ^ 0, in case q > 0. 



Since a > —1/2, it follows that each entry of M has an integrable singularity at the origin. Because 
M(0 = 0(1/^^) as C — *■ OO) and M is analytic in the upper half plane, it then follows by Cauchy's 
theorem that f^M+{C)dC = 0, and hence by 



/ H+{C)H-{CrdC = 0. 

JR 

Adding this equation to its Hermitian conjugate, we find 



(2.22) 



Using dUni), and the fact that (e*(l^'+<)'^«)* = e-*(t^'+'^^)'^3 for C, s G M, (here we use the 

fact that s is real!), we obtain from ()2.22|) . 



F_(C)(o l]H_iCrdC = 2l \iH^),2iO\' + \iH-)22iO\' 



dC. 
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This implies that the second column of H_ is identically zero. The jump relations H2.17() and H2.18() 
of H then imply that the first column of is identically zero as well. 

To show that the second column of and the first column of H_ are also identically zero, we 
use an idea of Deift et al. ^1 Proof of Theorem 5.3, Step 3]. Since the second column of H- is 
identically zero, the jump relations 1)2. 17() and (|2.18() for H yield for j = 1, 2, 

{Hj2)+ (C) = -e^g'^(^)-*-e-2*(t^'+<) iH,i)_ {(), for C G M \ {0}. 
Thus if we define for j = 1,2, 

p,.(C), iflmOO. 
[H.dO. ifImC<0, 

then both hi and /i2 satisfy the following RH problem for a scalar function h. 

RH problem for h: 

(a) h is analytic on C \ M and satisfies the following jump relation 

h+{C) = _e«g°K)™e-2^(tf'+<)/i_(C), for C G M \ {0}, 



(b) h{0 = O{l/0 as C ^ oo 

(c) h{C) 



fodCr), as C ^ 0, in case a < 0, 
[Oder"), as C ^ 0, in case a > 0. 



Take Co with Im^o < — 1 and define 

^^/i(C), ifImC>0, 



HO 



^ (-e™e-2^(tf'+<)) h{C), if -1 < ImC < 0, 



(2.24) 



where we use principal branches of the powers, so that is defined with a branch cut along 
the negative real axis. Then it is easy to check that h is analytic in lm(^ > —1, continuous and 
uniformly bounded in Im C > — 1 , and 

h{() = 0(6"^'^*^^'^), as C ^ oo on the horizontal line ImC = — 1- 



By Carlson's theorem, see e.g. |31) . this implies that /i = 0, so that /i = 0, as well. This in turn 
implies that hi = and h2 = 0, so that H = 0. Then also ^' = and the proposition is proven. □ 

As noted before, Proposition 12.51 has the following consequence. 

Corollary 2.6 The RH problem for \I'o,, see Section 2.1, has a unique solution for every s G M 
and a > —1/2. 
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2.3 Proof of Theorem fTTl 



Theorem 11.11 follows from the connection of the RH problem for of Section 2.1 with the RH 
problem associated with the general Painleve II equation (|1.12j) as first described by Flaschka and 
Newell HH! Section 3D]. 



Proof of Theorem II. XL Consider the matrix differential equation 



dc 



(2.25) 



where A is as in (|1.16j) and s, q, and r are constants. For every k = 0,1, ... ,5, there is a unique 
solution ^-fc of such that ^-^(0 = (I+0(l/C))e-*(5<'+0^3 as C ^ oo in the sector (2A:-l)f < 

argC < (2A; + 1)|. The function 

^(C) = ^fc(C), for (2A:-l)f <argC< (2A: + l)f, (2.26) 

is then defined on C \ (S U iM) and satisfies the following conditions. 

(a) ^' is analytic in C \ (S U iM). 

(b) There exist constants si, 82,83 £ C (Stokes multipliers) satisfying 

•si + S2 + S3 + 818283 = — 2zsin7ra (2.27) 
such that the following jump conditions hold, where all rays are oriented to infinity, 




on Fi 



on 



on Fo 




on Fq 



on 



on F4 



(c) ^(C) = (/ + 0(l/C))e-*(t^'+"^)"^ as C 



00. 



The Stokes multipliers si, S2, S3 depend on s, q and r. However, if q = g(s) satisfies the second 
Painleve equation q" = sq + 2q^ — a, and if r = q'{8), then the Stokes multipliers are constant. In 
this way there is a one-to-one correspondence between solutions of the Painleve II equation and 
Stokes multipliers 81,82,83 satisfying ()2.27|) . This also means that there exists a solution of the 
above RH problem which is built out of solutions of ()2.25|) if and only if s is not a pole of the 
Painleve II function that corresponds to the Stokes multipliers 81,82,83. The Painleve II function 
itself may then be recovered from the RH problem by the formula ^Hl 

q{8) = lim 2iC^i2(C)e-*(t^'+'^\ 

with ^'12 the (1, 2)-entry of ^ . In particular, condition (c) of the RH problem can be strengthened 
to 



*(C) 



/ + 



1 
2iC 



u(s) 
-q{8) 



-u{8) 



+ 0{l/e)]e 



as C ^ 00, 



(2.28) 
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where u = (g')^ ~ ■^q'^ ~ ^ 2aq. 

The RH problem for ^'^ in Section 2.1 corresponds to 

= e"™, S2 = 0, S3 = -e™. (2.29) 

These Stokes multiphers are very special in two respects |^[^. First, since S2 = 0, the corre- 
sponding solution of the Painleve II equation decays as s ^ +oo, i.e., 

q(s) ~ — , as s — > +00. (2.30) 

s 

Secondly, since sis^ = —1 the Painleve II solution increases as s — > — oo, i.e., 

q{s) ~ —- , as s — > — oo. (2-31) 

where the choice si = e~'^*", S3 = — e'^*" corresponds to the +sign, while the interchange of si 
and S3 corresponds to the — sign in (|2.3H) . Thus the special choice (|2.29|) corresponds to Qa, the 
Hastings-McLeod solution of the general Painleve II equation, see H1.13() and 1)1.14(1 . 

Then as a consequence of the fact that the RH problem for stated in Section 2.1 is solvable 
for every real s by Corollary 12.61 we conclude that Qa has no poles on the real line, which proves 
Theorem 11.11 □ 



Remark 2.7 Its and Kitaev [2l| use a slightly modified, but equivalent, version of the RH problem 
for The solutions are connected by the transformation 

^a^e'^'^'^^ae-^"^, (2.32) 
which results in a transformaton of the Stokes multipliers Sj (— l)-'zsj. 

For later use, we record the following corollary. 

Corollary 2.8 For every fixed sq € M, there exists an open neighborhood U of sq such that the RH 
problem for is solvable for every s G U. 

Proof. Since Qa is meromorphic in C, there is an open neighborhood of sq without poles. This 
implies jl6j that the RH problem for is solvable for every s in that open neighborhood of sq, as 
well. □ 

Remark 2.9 The function ^aiC'i s) is analytic as a function of both ^ € C \ E and s € C\Va, 
where Va denotes the set of poles of qa, see ^HI- As a consequence, one can check that (|2.5() . (|2.6() 
and p.7|) hold uniformly for s in compact subsets of C \ . 

Remark 2.10 The functions ^a,i and ^a,2 defined by (jl.l5|) and p. 18(1 are connected with as 
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follows. Define 



1 0' 



1 

.0 1 ) 

.0 1 , 



1 0' 

-Tvia 

1 0^ 



for C G 5i, 
for C G ^2, 
for C G ^3, 

for C G 5'4, ReC > 0, 

for C G 54, ReC < 0. 



(2.33) 



Then it follows from the RH problem for that <I>q, is analytic on C \ (—ioo, 0]. Morever, we also 
see from (|1.15p and ()1.18p that 



a,2 



(2.34) 



where * denotes an unspecified unimportant entry. It also follows that ^a,i and ^a,2 have analytic 
continuations to C \ (— ioo,0]. 

Remark 2.11 We show that the kern K^^^^'°'{u,v; s) is real. This will follow from the identity 



e^™"snW$^_2(u;s) = e^™"S'^(")«>„,i(u;s), 
since obviously ()2.35|) implies that 



for u G M \ {0} and s G 



(2.35) 



K''"''"{u,v;s) = — -Im f e5™('^s>^(")-^s°(^))$« i(n; s)$a l(^'; s)) • 

7r[u — V) \ ' ' / 



The identity (|2.35j) will follow from the RH problem. It is easy to check that ai^aiC] s)ai, with 
ai = ill), also satisfies the RH conditions for ^'q, . Because of the uniqueness of the solution of 
the RH problem, this implies 



^'„(C;s) = cTi^'„(C;s)fTi. 
For C G S4, the equality of the (2, 1) entries of yields by (ITHHl) and dOl) 



e"*"$a,2(C;s) = ^a,i(C;s), for C G 54, ReC > 0, 



and 



e-™$a,2(C;s) = ^a,i(C;s), for C G 54, ReC < 0. 



(2.36) 
(2.37) 
(2.38) 



Since both sides of ()2.37p are analytic in the right half-plane we find the identity 1)2. 35(1 for n > 0, 
and similarly since both sides of (|2.38|) are analytic in the left half-plane, we obtain (|2.35j) for u < 0. 
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3 Steepest descent analysis of RH problem 



In this section we write the kernel Kn,N in terms of the solution Y of the RH problem for orthogonal 
polynomials (due to Fokas, Its and Kitaev |17j l and apply the Deift/Zhou steepest descent method 
jl5j to the RH problem for Y to get the asymptotics for Y . These asymptotics will be used in the 
next sections to prove Theorems 11.21 and 11.71 

We will restrict ourselves to the one-interval case, which means that ipv is supported on one 
interval, although the RH analysis can be done in general. We comment below in Remark 13. II fsee 
the end of this section) on the modifications that have to be made in the multi-interval case. 

As in Theorems 11.21 and 1 1 . 71 we also assume that besides there are no other singular points. 

3.1 RH problem for orthogonal polynomials 

The starting point is the RH problem that characterizes the orthogonal polynomials associated 
with the weight \x\'^°'e~^'^^^\ The 2x2 matrix-valued function Y = Y^^n satisfies the following 
conditions. 



RH problem for Y 

(a) y : C \ M ^ C2^2 ig analytic. 

'l |xP°e"^^('') 



(b) y+(x) = y_(x) 







1 



for X E M. 



(c) Y{z) = {I+0{l/z)) 







as z ^ oo. 



" 

(d) y has the following behavior near the origin 



Yiz) 




as z ^ 0, if a < 0, 



as z — > 0, if a > 0. 



(3.1) 



Here we have oriented the real axis from the left to the right and y+(x) {Y-{x)) in part (b) denotes 
the limit as we approach x E M from the upper (lower) half-plane. This RH problem possesses a 
unique solution given by \^t\ (see HTjUHl for the condition (d)). 



Y{z) 



( 



\ 



1 



1 



-2mKn-l^NPn-l,N{z) -Kn-l,N 



— dy 



2lTiKn,N Jr y - Z 

p„_i,^(y)|y|2°e-^^(^) 



y - z 



dy 



(3.2) 



for z E C\M, where Pn,N{z) = «^n,Af-2" + ' • • , is the n-th degree orthonormal polynomial with respect 
to the weight \x\^°'e~^'^^^\ and Kn,iV is the leading coefficient of Pn,N- 
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The correlation kernel K^^n can be expressed in terms of the solution of this RH problem. 
Indeed, using the Christoffel-Darboux formula for orthogonal polynomials, we get from ()1.4I) . H3.2() . 
and the fact that dety = 1, 

J^n,N\-i'^ y — c ^ \y\ e ^ 

fin,N X-y 



I |ag-|7Vy(a;)| |ag-l7VV{y) 



The asymptotics of K^^n follows from a steepest descent analysis of the RH problem for Y , 
see (71 El El UHl Uni ES] • The Deift/Zhou steepest descent analysis consists of a series of explicit 
transformations Y\-^T^S^R so that it leads to a RH problem for R which is normalized at 
infinity and which has jumps uniformly close to the identity matrix I. Then R itself is uniformly 
close to /. By going back in the series of transformations we then have the asymptotics for Y from 
which the asymptotics of Kn,N in different scaling regimes can be deduced. 

The main issue of the present situation is the construction of a local parametrix near with 
the aid of the RH problem for introduced in Section (21 For the case a = this was done in |7] 
and we use the ideas introduced in that paper. 

Throughout the rest of the paper we use the notation 

t = n/N, and ^i = \^- (3.4) 
3.2 First transformation Y ^ T 

In the first transformation we normalize the RH problem at infinity. The standard approach would 
be to use the equilibrium measure in the external field Vt, see [HI I32j- This is the probability 
measure that minimizes 

IvAlj) = j j log 1^ \ y^ dfi{x)dfj.{y) + j Vt{x)dfi{x) 

among all Borel probability measures fi on M. The minimizer for t = 1 has density ipv which by 
assumption vanishes at the origin. For t < 1, the origin is outside of the support and for t slightly 
less than 1, there is a gap in the support around 0. An annoying consequence is that the equality 
in the variational conditions is not valid near the origin. Therefore, a modified measure vt was 
introduced in to overcome this problem. 

Here, we follow O Section 3]. We take a small (5o > so that tpv{x) > for x £ [—Sq, 6o] \ {0}, 
and we consider the problem to minimize IvA^^) among all signed measures v = — where 

are nonnegative measures such that f du = 1 and supp(z^~) C [— (5o,(5o]. There is a unique 
minimizer which we denote by Vf This signed measure is absolutely continuous with density ^pt 
and its support St = [at , bt] is an interval if t is sufficiently close to 1 . The following variational 
conditions are satisfied: there exists a constant G M such that 

2 j log \x - y\My)dy - Vt{x) +it = 0, for x € [at, h], (3.5) 

2 j log [x - y[^l)t{y)dy - Vt{x) + ft < 0, for x G M. (3.6) 
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In addition, it was shown in [7j that for t sufficiently close to 1, 
Mx) = for X G [at, h], 

where 

■V\z)Y 1 fV'{z)-V'{y) 



Qt{z) 



2t 



z-y 



-ipt{y)dy. 



(3.7) 



(3.8) 



For t > 1, we take the square root in (|3.7|) which is positive for x = 0, while for t < 1 we take the 
square root which is negative for x = 0. 

For the first transformation, we introduce the following '(7- function' associated with vt, 

gt{z) = j \og{z - y)dut{y) = j log(z - y)My)dy, for z G C \ M, (3.9) 

where we take the branch cut of the logarithm along the negative real axis. We define 

T{z) = el"^*'^3y(^)g-n9tW<x3g-|n£ta3^ for z G C \ R. (3.10) 
We also use the functions 

(3.11) 



Jbt 



J at 



/^ds, 



(3.12) 



where the path of integration does not cross the real axis. The relations that exist between gt, (ft 
and (ft are described in |7|l Section 5.2]. Using these, we find that T is the unique solution of the 
following RH problem. 

RH problem for T: 

(a) T : C \ M ^ C2x2 is analytic. 

(b) r+(x) = T_{x)vt{x) for x G M, with 

1 2a \ 



vt{x) = < 











1 



'1 |3;|2ag-2n<^t(a::)^ 







1 



for X G iat,bt), 
for X G {bt, 00), 
for X G (—00, at). 



(c) T{z) = 1 + 0{l/z), as z ^ 00. 

(d) T has the same behavior as Y near the origin, given by (|3.1|) . 
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Figure 2: The lens shaped contour Eg going through the origin. 



3.3 Second transformation S t-^ T 



In this subsection, we open the lens as in Figure [2 The opening of the lens is based on the 
factorization of the jump matrix vt for x £ (at, bt), which is 



vt{x) 







\x 



2a 



\^\-2a^2rupt~{x) ]^ 



2a 



-2a 







\^\-2a^2rnpt,+ {x) 



(3.13) 



We deform the RH problem for T into a RH problem for S by opening a lens around [at, ht] going 
through the origin, as shown in Figure |2I The precise form of the lens is not yet specified but for 
now we choose the lens to be contained in the region of analyticity of V and we can do it in such a 
way that for any given 5 > 0, there exists 7 > so that, for every t sufficiently close to 1, we have 
that 



Yieipt{z) < -7, 



(3.14) 



for z on the upper and lower lips of the lens with the exception of (5-neighborhoods of 0, a, and b. 
See also Section 5.3]. 

Let LO be the analytic continuation of x 1— > \x\'^°^ to C \ («M), i.e.. 



uj{z) 



ifRe2>0, 
(-2)2", ifRe2<0. 



(3.15) 



The second transformation is then defined by 
't(z), 

T{z) 



S{z) 



Tiz) 



1 0' 

^-(j(z)-ie2"^*(^) 1^ 

' ' 



for z outside the lens. 

for z in the upper parts of the lens, 

for z in the lower parts of the lens. 



(3.16) 



Then S is the unique solution of the following RH problem posed on the contour which is the 
union of M with the upper and lower lips of the lens. 
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RH problem for S: 

(a) 5 : C \ S5 ^ C2X2 is analytic. 

(b) S+ = S-Vs on 'Es, where 
■ ' 1 0' 

|2a\ 



vsiz) = < 



\z\ 
,-kr'° ) ' 

^1 |^|2ag-2n^t{2)^ 

1 



for z G Ss \ K, 
for 2: G iat,bt), 
for z G (6t, 00), 
for z G (—00, at). 



(c) S'(z) = / + 0(1/2), as z ^ cx). 

(d) has the following behavior near the origin. If a < 0, 



1 z 



2a 



S{z)=0[^ ! 2a), asz^O,zGC\Ss, 



and if a > 0, 



S{z) = < 



O 



1 1 



as z — > from outside the lens, 



/|z|-2" i\ 

O \ „ , as z — > from inside the lens. 

\lz|-2° ir 



(3.17) 



(3.18) 



3.4 Parametrix away from special points 

On the lips of the lens and on (—00, at) U (6t,oo), the jump matrix for S is close to the identity 
matrix if n is large and t is close to 1. This follows from the inequality 1)3. 14() and the fact that 
ftix) > for X > bt and ^ti^) > for x < at- Ignoring these jumps we are led to the following RH 
problem. 

RH problem for 

(a) P(°°) ■.C\[at, bt] C2x2 is analytic. 



(b) P|°°)(x) =p!°")(x) 

(c) P(°°)(z) =/ + 0(1/z), 



\x 



2a 







as z — > 00. 



for X G {at,bt) \ {0}. 
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Note that depends on n and N through the parameter t. As in j25l 1271 OH] we construct 

p(°o) in terms of the Szego function D associated with on {at,bt). This is an analytic function 
in C\[at, bt] satisfying D^{x)D^{x) = for x £ {at, bt)\{0} and which does not vanish anywhere 
in C \ [at, bt]- It is easy to check that D is given by 

I)(z) = z>f ^\~°*~^* ') °, for zeC\[at,bt], (3.19) 
V bt - at ) 

where ^{_z) = z + {z — + lY^"^ is the conformal map from C \ [—1, 1] onto the exterior of the 

unit circle. Since = 2z + 0{1/ z) as z ^ oo we have 



\\mD{z)=\- =D^. 

z^oo yij^-atj 

Now the transformed matrix valued function 

= i)-<^3p{oo)^a3 (320) 

satisfies conditions (a) and (c) of the RH problem and it has the jump matrix o) on {at,bt). 
The construction of has been done in |51 El leads us to the solution of the RH 

problem for 



-1 



' f^{z)+f^{z)-^ I3{z)-^{z) 

P'^°^\z) = D^\ ^ , 1 ior zeC\\at,bt], (3.21) 

^ ^ °° ' f3{z)-l3{z)-^ I3{z)+I3{z)-^ I ^ ^ ' \ I t) tj) \ J 



-2i 2 



where 



Piz) = l' iovz€C\[at,bt]. (3.22) 

[z - 04)^4 



3.5 Parametrix near endpoints 

The jump matrices of S and P^°°^ are not uniformly close to each other near the origin and near 
the endpoints of [at,bt]. We surround oi and 61 (the endpoints of Sy) with small disks Us{a) and 
Us(b) of radius 5. For t sufficiently close to 1, the endpoints at and bt are in these disks, and then 
local parametrices and can be constructed with Airy functions as in [51 1131 ITl] . 



3.6 Parametrix near the origin 

Near the origin a local parametrix will be constructed with the aid of the RH problem for ^q, of 
Section 121 Let Us be a small disk with center at and radius 6 > 0. We seek a 2 x 2 matrix valued 
function P in Us with the same jumps as S, with the same behavior as S near the origin, and which 
matches with p(°°) on the boundary dUs of the disk. We thus seek a 2 x 2 matrix valued function 
that satisfies the following RH problem. 

RH problem for P: 

(a) P is defined and analytic in Us' \ S5 for some 6' > 6. 
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(b) On T,s n Us, P satisfies the jump relations 

P+{z) = P.{z) f ,,Jmz) ?1 . for z G (Ss \ M) n Us, (3.23) 



P+(x)=P_(x)(_|^°_2„ '"^l^'"), forxGM,5)\{0}. (3.24) 

(c) P satisfies the matching condition 

P{z) = {I + 0{n-'^/^))P^^\zy, (3.25) 
as n, ^ cxD such that n'^/'^{n/N — 1) — > L, uniformly for z S dUs \ S5. 

(d) P has the same behavior near the origin as S, given by (|3.17() and 1)3.18(1 . 

In order to solve the RH problem for P we work as follows. First, we seek P such that it satisfies 
conditions (a), (b), and (d). To do this, we transform (in the first step) the RH problem into a RH 
problem for P with constant jump matrices. In the second step we solve the RH problem for P 
explicitly by using the RH problem for ^'q. In the final step we take also the matching condition 
(c) into account. 

Step 1: Transformation to constant jump matrices 

In the first step we transform the RH problem for P into a RH problem for P with constant jump 
matrices. We seek P in the form 

P{z) = S(z)P(z)e"^*(^)'"^e5™'^3^-°'^3 ^ if ^ > 0, (3.26) 

P{z) = ^(z)P(z)e~"^*(^)'"»e^™'"3 (^J if Imz < 0, (3.27) 

where the invertible matrix valued function E = En,N (we suppress notation of the indices) is 
analytic in Us' and where the branch cut of z" is chosen along the negative real axis. 

Using (|3.24j) . ()3.26|) and (|3.27|) . and keeping track of the branches of we can easily check that 



J' 



P has no jumps on {—S, 5) \ {0}. What remains are jumps on the contour (S5 \M) H C/^ = Uj=i ^ 
which is shown in Figure El We have reversed the orientation of S2 and S3 towards infinity, so that 
now the orientation of the Sj's corresponds to the orientation of the Fj's in Figure ^ The contour 
divides Us into four regions I , II , III and IV , also shown in Figure IHl 

We will now determine the jump relations for P. By ((3.15|) . (|3.23|1 . and (|3.26|1 . P should have 
the following jump matrix on Si, 



1 o\ _ / 1 



^aas g- iTriacra ^-nift (2)0-3 



(3.28) 



For z G S2 we have, because of the reversal of the orientation, an extra minus sign in the (2, l)-entry 
of the jump matrix. The result is 

P-{z)-'P+{z) = (^_^^^^J^,,^_^,^ ?) = (-e™ l)' ^^'^^^ 
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Figure 3: Contour and jumps for the RH problem for P. 



where the last equality follows from the fact that u){z)^^z'^°' = e^'^*" by ()3.15l) . since Rez < in this 
case. Using equations (|3.24|) and (|.S.27|) . the jump matrices for P on S3 and S4 can be determined 
similarly. The result is that 



P-iz) 



1 e-""' 
1 



1 



for z G S3, 
for z £ S4. 



(3.30) 



We arrive at the following RH problem for P. If it is satisfied by P then P defined by (|3.2(i)l - 
H3.27() satisfies the parts (a), (b), and (d) of the RH problem for P. 

RH problem for P: 

(a) P is defined and analytic in Us> \ [jj Sj for some 5' > 6. 

(b) P satisfies the following jump relations 



P+iz) = { 



P~{z) 

P-{z) 
P.{z) 

P-{z) 



1 

,0 1 , 



for 2; G Si, 
for z G S2, 
for z G S3, 
for z G S4. 



(3.31) 
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(c) P has the fohowing behavior near the origin. If a < 0, 



P{z) = O 



and if Q > 0, 



as z — > 0, 



(3.32) 



P{z) = I 



O 



O 



o 



z|" \z\ 

\z\ 



j , as z ^ 0, z G I U III , 
, as z ^ 0, z G II , 
, as z ^ 0, z G IV . 



(3.33) 



Note that if P has the behavior near the origin as described in part (c) of the RH problem, then P 
defined by ()3.26|) and (|3.27|) has the same behavior near the origin as S, as required by part (d) of 
the RH problem for P. 

Step 2: Construction of P 

Observe that the jump matrices and the behavior near the origin of the RH problem for P corre- 
spond exactly to the jump matrices and the behavior near the origin of the RH problem for ^ a. 
We use the solution of the latter RH problem to solve the RH problem for P. 
We seek P in the form 

P{z) = {n^'^f{z)-n^'\st{z)) , (3.34) 

where / and st are analytic functions on Us which are real on {—5,8), and sj is such that 

r?/'^st{z) eC\Va, for z G Us, (3.35) 

where Va is the set of poles of Qa- In addition, / is a conformal map from Us onto a convex 
neighborhood f{Us) of such that /(O) = and /'(O) > 0. Depending on / we open the lens 
around [at, bt] such that that /(Sj) = Fj for i = 1,2, 3, 4, where the Fj's are the jump contours for 
the RH problem for \I'q, see Figure ^ Recall that the lens was not fully specified and we still have 
the freedom to make this choice. 

It remains to determine / and st so that the matching condition for P is also satisfied. Here 
we again follow j7|. As in [3 Section 5.6] we take 



1/3 



vrV'y(O) 



1/3 



+ O (z^) , as 



0, 



and 



st{z)f{z)= 1^ ((-Q,(y))l/2_(_g^(y))l/2j^y_ 



(3.36) 



(3.37) 
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Then / is analytic with /(O) = and /'(O) > 0, it does not depend on t, and it is a conformal 
mapping on Us provided 6 is smah enough. Since the right-hand side of ()3.37|) is analytic and 
vanishes for z = 0, we can divide by f{z) and obtain an analytic function st- From [71 (5.29)], we 
get that there exists a constant K > such that 

\stiz) -TTC^^^ (t- l)wsy (0)1 < K{t-l)\z\ +o{t- I) ast^l, (3.38) 

uniformly for z in a neighborhood of 0. Now assume that |n^/^(t — 1)| < M and n large enough. 
Then it easily follows from ()3.38|) and the fact that qa has no real poles, that there exists a 6 > 0, 
depending only on M, such that 

\lmri^^^ st{z)\ < min{|Ims| | s is a pole of Qa} for \z\ < 5. (3.39) 

Then (|3.35j) holds and (|3.34|) is well-defined and analytic since ^aiCi is jointly analytic in its two 
arguments, see Remark 12.91 

It follows from (lOfill and that 



—I 



.J{zf + st{z)f{z) 



99t,+ (0) -ipt{z) 

(/?t,+ (0) +V9t(z) 



if Im z > 0, 
if Im z < 0, 



(3.40) 



see also Section 5.6]. Hence by (|2.5|) . which by Eemark 12.91 holds uniformly for s in compact 
subsets of C \ "Pq , we have 



P{z) = M/, ( n^'^f{z)-n"'\st{z) ) I + 0(l/n^/^) 



2/3, 



l/3^ 



„n</3t, +(0)0-3 



X < 



g-nv^t(z)(73 iflniz>0, 



if Im z < 0, 



(3.41) 



as n, — > 00, uniformly for z G dUs- 



Step 3: Matching condition 

In the final step we determine E such that the matching condition (c) of the RH problem for P is 
satisfied. By (jS^Hl), (E2ZI), and (HT:!!]) we have for z e dUs, 



P{z) 



if Im z > 0, 



l™<73 f 1 1 ^-aa3 if Iniz < 0, 



1 



as n, ^ CO. This has to match the outside parametrix p(°°\ so that we are led to the following 
definition for the prefactor E{z), for z £ Us, 



' (2) Z°'^3 g- ^TT-iaaa ^-nipt,+ (0)0-3 ^ 



E(z) 



if Im 2; > 0, 



1 



p(oo)^^^^a<73 ( M g-i7riao3g-n¥,i,+(0)fr3^ if Im Z < 0. 



(3.42) 



One can check as in [291 134j that E is invertible and analytic in a full neighborhood of Us- In 
addition we have the matching condition 1)3. 25() . This completes the construction of the parametrix 
near the origin. 



26 



Figure 4: The contour S/j after the third and final transformation. 



3.7 Third transformation: S t-^ R 

Having the parametrices P^°°\ p("-t)^ and P, we now define 

'S{z)P-\z), iovzeUs, 

5(z)(p(^*))^'(z), fovzeUsib), 
^S{z) (PM)"' (z), for z G C \ (Us U Usia) U Usib) U 



(3.43) 



Then R has only jumps on the reduced system of contours S/j shown in Figure |1J and R satisfies 
the following RH problem, cf. J_^. The circles around 0, at and bt are oriented counterclockwise. 



RH problem for R: 

(a) P : C \ ^ C2^2 ig analytic. 

(b) R+{z) = R-{z)vfi{z) for z G S/j, with 
>(oo)(p(aO)-i^ ona[/5(a), 

p(oo)(p(fet))-i, ondUsib), 

p(^)p-\ on 

^p(°o)^;_5(p{°o))-i, on the rest of T,r. 

(c) R{z) = 1 + 0(1/2), as 2: ^ 00, 

(d) R remains bounded near the intersection points of E^j. 



VR= \ 



(3.44) 



Now we let n, ^ cxD such that \n'^^'^(n/N — 1)| < M, so that 6 does not depend on n. Then 
it follows from the construction of the parametrices that 

'l + 0{l/n), on dUs (a) UdUsib), 
VR= h + 0{n-^/''), on dUs, (3.45) 
/ + ©(e-T"), on the rest of Sr, 

where 7 > is some fixed constant. All O-terms hold uniformly on their respective contours. 
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For large n, the jump matrix of R is close to the identity matrix, both in L°° and in L^-sense 
on Tiji. Then arguments as in 1131 114j (which are based on estimates on Cauchy operators as 
well as on contour deformations), guarantee that 

Riz) = I + 0{n-^^^), uniformly for z e C \ Sij, (3.46) 

as n, N ^ oo such that \n'^/^{n/N - 1)| < M. 

This completes the steepest descent analysis. Following the effect of the transformation on the 
correlation kernel Kn^N and using ()3.46|) we will prove the main Theorem 11.21 This will be done 
in the next section. For the proof of Theorem 11.71 we need to expand vji{z) in (|3.45j) up to order 
n~^/^, from which it follows that 

R(z) =1+ ^^-l^ + 0(n-2/3), uniformly for zeC\^R, 

with an explicitly computable R^^\z). The asymptotic behavior of the recurrence coefficients is 
expressed in terms of i?^^) and this leads to the proof of Theorem II. 71 This will be done in Section 

ini 

Remark 3.1 The steepest descent analysis was done under the assumption that supp('i/'v') consists 
of one interval. In the multi-interval case, the construction of the outside parametrix p(°°) is more 
complicated, since it uses 0-functions as in ^1 Lemma 4.3] and the Szego function for multiple 
intervals as in [29' Section 4]. With these modifications the asymptotic analysis can be carried 
through in the multi-interval case without any additional difficulty. 

4 Proof of Theorem 11.21 

As in the statement of Theorem 11.21 we assume that n, N ^ oo with in?^^{t — 1) ^ L, where 
t = n/N . Let M > \L\ and take n sufficiently large so that |n^/^(t — 1)| < M. Let 5 > be such 
that H3.39() holds. We start by writing the kernel Kn,N explicitly in terms of the matrix valued 
function defined in (|2.33|) . For notational convenience we introduce 

B{z) = R{z)E{z), (4.1) 

where E and R are given by (|3.42() and ()3.43|) . respectively. 

Proposition 4.1 Let x,y G (—(5,5) \ {0}. Then 

= 27ri(x-j/) """"^'"'^''^^''°^'^^ (° ^) ("'/'/(y); n'^^stiy)) 

X B-\y)B{x)'^^ {n'/'f{xy,n'/ht{x)) Q , (4.2) 

where "is given by (|2.33|) . 

Proof. From (|3.3|) . (|3.1()|) . and the fact that NV = nVt, the kernel Kn^N can be written as 
A'„,K(i,!/) = |i|°e3"<2''"<'^>-'''<'"'+*'|j/rel»<^»''+W-''''»'+''> 

"2^,^' l)T-fa)T,w(j), 
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Using the relation 

2gt^+ -Vt+it = -2(pt,+ on [at, h], 
see (Tj, and (|3.16|) to express T in terms of S, we find for x and y in {at, bt) \ {0}, 

27rz(x — y) 



X 5+(x)|xr"3e-"^*-+(^)'^3 (^1^ . (4.3) 



We further simphfy this expression by writing S in terms of R and the parametrix P near the 
origin. Consider the case that x £ (0,(5). Then, since S+{x) = R{x)P^{x) by (|3.43j) . we have by 

S+{x) = S(x)P(x)e^™'"3e"'^*'+(^)'^3|2.|-a<x3^ ^ ^ (0,5), (4.4) 

where B is given by (gUJ). By (|OH) . and (fT^ we then find for x G (0, 5), 

5+(x)|xr"^e-"^*-+(^)"3 

= i?(x)$„ (nV3/(^);n2/3,,(x)) J) e^™'^3 (^ij 

= e^™^g°(^)s(x)$„ (n^/3/(x);n2/3st(x)) Q . (4.5) 

A similar calculation shows that (|4.5|) also holds for x G {—6,0). Similarly, we have 
(-1 1) |?/r"'^3en^t,+(?/)^3 5-i(y) 

= el^^s'^^^) (0 1) <P-' {n'/'f{yy, n'^ht{y)) B-\y), (4.6) 

for y S {—5,5) \ {0}. Inserting ()4.5p and 1)4. 6(1 into 1)4. 3() . we arrive at (|4.2p . which proves the 
proposition. □ 



Proof of Theorem 11.21 Let u,v GM. \ {0}, and put n„ = u/{cn^/^) and z;„ = t;/(cn^/3) with c 
given by THWi . Note that, by (HOHl) . 

lim n^/V(^in) = lim n^/V('^n) = (4.7) 
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Furthermore, by (E3HI), dH), and THM . 

\n'^/ht{z) -s\< Kn^/^{t - l)\z\ + n^/'^o{t - 1) + \n'^/^{t - 1) - L\Trc^/^wsy{0) 
uniformly for z in a neighborhood of 0. Then it easily follows that, since n 

lim n^/^st(n„) = lim n^^^stivn) = s. (4.8) 

n,N—*oo n,N—*oo 

Now, similar as in [2^], we use the fact that the entries of B are analytic and uniformly bounded 
in Us, to obtain 

lim B-\vn)B{un) = 1. (4.9) 
Inserting (|I?7|) . and into (g^, we find that 

lim —^Kn,N{Un,Vn) 
n,Af— >oo cn ' 



27rz(ti — v) 



,|™(sgnH+sgnW) ^0 1) s)<^a{u; s) 



2-Ki{u — v) 

This completes the proof of Theorem II. 21 □ 

5 Proof of Theorem 11.71 

In this section we will determine the asymptotic behavior of the recurrence coefficients an^N and 
bn^N sls n, N ^ oo such that \n'^^^(n/N — 1)| < M for some M > 0. As in Theorem 1 1.71 we assume 
that Sy = [a, b] is an interval, and that there are no other singular points besides 0. Then it follows 
that supp('(/'t) consists of one interval [at, bt] if t is sufficiently close to 1. In addition we have that 
the endpoints at and bt are real analytic functions in t, see |26| Theorem 1.3], so that 

at = a + 0(n-2/3), bt = b + 0{n'^/^), (5.1) 

since t = n/N = 1 + ©(n'^/S). 

We make use of the following result, see for example |H1E|- Let Y be the unique solution of the 
RH problem for Y. There exist 2x2 constant (independent of z, but depending on n, A^) matrices 
Yi , Yz such that 



and 



(i^2) 



an,N = \/(yi)l2(yi)21, bn,N = (^l)ll + 77^- (5.2) 

We need to determine the constant matrices Yi and Y2. For large |z| we have by 1)3. 1U() . ()3.16|1 and 

that 

Y{z) =e-^"^""^ii(z)p(°°)(z)e"f*(^)'"3ei"^"^3_ (5,3) 

So in order to compute Yi and Y2 we need the asymptotic behavior of P^°°\z), e"'^*^^)'^^ and R{z) 
as z ^ cxD. 
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Asymptotic behavior of P^^\z) as z ^ oo: 

From ()3.19|) and ()3.22|) it is straightforward to determine the asymptotic behavior of the scalar 
functions D[z) and f3{z) as z — > oo. Indeed, as z ^ oo, 

2 2i 

I3{z)-f3{z)-^ I3{z)+p{z)-^ 
-2i 2 



—i\ 1 i ,,0 2n 



* 1\ 1 



-1 * } z"^ 



+ 0{l/z- 



and 



D{z)- 



a 



I--{ht + at) 



1 \ 1 /* 0\ 1 



-II z^ 



„ - + 0{l/z^) 
* J z^ 



where * denotes an unspecified unimportant entry. Inserting these equations into (|3.21|1 and using 
(|5.H) gives us the asymptotic behavior of at infinity. 



p(oo) p{oo) 



(5.4) 



with 



(oo) 



-^{bt + at) lih- at) 



oo 



a 



2V- + «) lib -ay 
1(6 -a) f(6 + a)^ 



Z)^-3+0(n-2/3), 



(5.5) 



and 



(~) 



oo 



(a-l)(62-a2) 



^ oo 



.(a + l)(62-a2)^ 



Asymptotic behavior of e"^*^^)'^^ as z — > oo: 
By ()3.9|) we have 

jngt{z)az / ^ " 



z' 



I + ^ + ^ + 0{l/z^) 



z z^ 



Z)-'^«+0(n-2/3). 



as z — > oo, 



(5.6) 



(5.7) 



with 



Gi = -n yiJt{y)dy 



1 

-1 



Go 



* 
* 



(5.8) 
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Asymptotic behavior of R{z) as z ^ oo: 

The computation of Ri and R2 is more involved. For z € dUs H C+, we have by H3.44() . H3.26() . 
^M), and ((321, 

vr{z) = P^^\z)p-\z) 

Using (fT^ and (HTini) . we then find 

vr{z) = I + — ^ + 0(n-2/3), (5.10) 



n 



where 



-1 

^(l)(^-) _ p(oo) ( ~) -"""a ^7riao-3 g-n(pf .+ (0)0-3 

2^/(2) 

\-qa{n'^/^St{z)) -Ua{n^/^st{z))) 

for z G OJ/i n C+. A similar calculation leads to an analogous formula for z G dUs H C_, which 
together with 1)5.11(1 shows that A^^^ has an extension to an analytic function in a punctured 
neighborhood of with a simple pole at 0. 

To calculate the residue at 0, we use (|3.19p together with the fact (/>+(x) = exp(z arccos x) for 
X G [-1, 1] to find 

D{z) ( . f h + at 

lim = exp I —la arccos 



so that by ^TT)^ 

lim p(°°)(^)2"'^3g-|™<x3 ^ D^p[°°\0)e''^'^"''', with 9t = arcsin (5.12) 



Also note that by (EHI), (EUl), (|T3!I|) and 

rh 

—(pt^+{0) = TTt il)t{x)dx = -irioJt. (5.13) 
Jo 

Now use Xrm . (Eini), and (|KTT|) in to find 

Res(A(i);0) = _-Li5^3p|~)(o)e*(^"'^*+"^*)'^3 



n„(n2/3s,(0)) g„(n2/3,,(0)) ^ 
-ga(n2/3st(0)) -u«(n2/3st(0))^ 
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Combining (lOTl) . and (glZI) we see that n'^/^st{0) = St,n as defined in (ITTTl) . From (HT^ . 

(jXTTl) . and (pr^ it follows that 



^ /3+(0)+/3+(0)-i /3+(0)-/3+(0)-n 

/3+(0)-|+(0)-^ /3+(0)+/3+(0)-i ' 
V 2j 2 / 



where /3+(0) = e^'^/^ (-6t/at)^/^ 

We insert this into (|5.14)) and after some straightforward calculations we find 

-Res(A«;0) = L»S(n^l+r2a2 + r3a3)^)^"^ (5.15) 

where the Pauli matrices are = ( ^ J ) > '^2 = ( ^ ^* ) , and 0-3 = ( q i'x ) ) and 

n = -7^ {ua{st,n) !^\ + qa{st,n)-;^^^==sm{2irnuJt + 206*4) ) 

= -7^ (ua{st,n):^^^ + ga(st,„)-^^^ sin(27rna;t + 2a^) ) + ©(n^^/^), (5.16) 
^^c \ 2v— ao 2\/—ab / 

r2 = ^^^^cos{27rnuJt + 2adt) 
2c 



qa{St. 



n. 



■ cos(27rnwi + 2ae) + ©(n'^/^), (5.17) 



2c 

'^s = 77- {qa{st,n) ^\ sin(27rna;t + 2a6'i) + nQ,(sj,„); ^* ^ 



= 7^ ( '7a(st,n)-^^ sin(2™i + 2a0) + n«(si^„)-l±^ ) + ©(n'^/^), (5.I8) 
2c V 2^J-ab 2-^ -ah) 

where we used 1)5. 1|) . 

From (f^lU)) it follows that 

i?(z) = /+^^^ + 0(n-2/3), (5.19) 

where R^^ = R^L^ + A^^) on dUs and i?(i)(z) ^ as 2; ^ 00. Since A^-^) is analytic with a simple 
pole at z = 0, we can find explicitly 

f-iRes(A«;0) + A«(z), for z G C/^, 
i?(^)(z) = < _ (5.20) 

[-iRes(A(i);0), for z G C \ t/^. 

As in jl4j the matrix valued function R has the following asymptotic behavior at infinity, 

R(z) = I+ — + ^ + 0{l/z^), asz^oo. (5.21) 
z z^ 

The compatibility with (|5.19|) and (|5.2Uj) yields that 

i?i = -Res(A(i); 0)n-^/3 ^ 0(^-2/3)^ = 0(n-2/3). (5.22) 

Now, we are ready to determine the asymptotics of the recurrence coefficients. 
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Proof of Theorem ITTTl Note that by B, (EIZl) and Km . 



and 



■k nit 0-3 



(5.23) 



(5.24) 



We start with the recurrence coefficient an,N- Inserting 1)5. 23() into (|5.2() and using (|5.5() and 
the facts that (Gi)i2 = (^1)21 = (by (|EH1)), and (iii)i2(iii)2i = ©(n-^/S) (by ^f^mh ). we obtain 



an,N 



1/2 



1/2 



From ()5.22|) and 1)5. 15() we then arrive at, 
h — a 



an,N 



-i/3 + 0(n-2/3) 



b-a _ qa{st,n) COs{2TTnuJt + 2a(9) „i/3 _2/3^ 

4 2c ^^'^ 



(5.25) 



Next, we consider the recurrence coefficient bn,N- Inserting (|5.23|) and (|5.24j) into (|5.2j) . and 
using the facts that (Gi)ii + (Gi)22 = (by (ESI), and (i?2)i2 = ©(n^^/S) (-^y ^^-^ obtain 



'Jn,N 



(P^°°))n + (i?i)ii + 



\ l-^l jl2 
^ / (^2 )l2 



{pr)i2+{RiPr)i2+o{n^h 

iPt^+Rl)l2 
^ (Pl)l2 ^ 



0(n-2/3) 



.IM J 12 



+ 



(Pi)n + ^^^(Pi)i2 + 0(n-^/^)). 
(^1 J 12 / 



From equations 1)5. 5(1 . 1)5. 6|) . ()5.22p . and ()5.15p . we then arrive at 

6n,JV = ^ + 2CRi 11 + 2i- Z^^' iJi 12 + O n-'/^' 

2 — a 



6 + a 



+ 2 (^rg + i^^(n - ir2)^ + ©(n^^/S). 



(5.26) 



Using 1)5. 16(1 . (|5.17p . and (|5.18p in (|5.26p we will see that the terms containing Mq, cancel against 
each other. What remains are the terms containing Qa- 



b + a Qaist,: 
bn,N = \ — 



b + a 2\/—ab . 

cos(27rnti;t + 206^) H — sm{2TTnuJt + 2aU) 



b — a 



b — a 



n 



"1/3 



+ 0(n-2/3). (5.27) 
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Since = sin 6 and = cos 9, we can combine the two terms within square brackets and the 

result is 

2 c 

Theorem [OH is proven by (|0^ and (|05|l . □ 
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